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We develop the microscopic theory of the extrinsic spin Hall conductivity of a two-dimensional 
electron gas, including skew-scattering, side-jump, and Coulomb interaction effects. We find that 
while the spin-Hall conductivity connected with the side-jump is independent of the strength of 
electron-electron interactions, the skew-scattering term is reduced by the spin-Coulomb drag, so 
the total spin current and the total spin-Hall conductivity are reduced for typical experimental 
mobilities. Further, we predict that in paramagnetic systems the spin-Coulomb drag reduces the 
spin accumulations in two different ways: (i) directly through the reduction of the skew-scattering 
contribution (ii) indirectly through the reduction of the spin diffusion length. Explicit expressions 
for the various contributions to the spin Hall conductivity are obtained using an exactly solvable 
model of the skew-scattering. 

PACS numbers: 
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I. INTRODUCTION 



There has recently been a strong revival of interest in the phenomenon of the spin Hall effect in the context of 
semiconductor spintronics"'^. Two different forms of this phenomenon have been identified: an extrinsic spin-Hall 
effect^-'^'^'^, which is driven by spin-orbit scattering with impurities, and an intrinsic onaSi^, which is due to spin- 
orbit effects in the band structure. In both forms, the phenomenon consists of the appearance of a transverse spin 
current - say a z-spin current in the y-direction - when the electron gas is driven by an electric field in the x-direction. 
The physical manifestations of this spin current are still an active subject of study and controversy.-^ However, it is 
now believed that the spin Hall effect should lead to transverse spin accumulation when the flow of the spin current is 
suppressed by an appropriate gradient of spin-dependent electrochemical potential. Recently, the spin Hall effect has 
been experimentally observed in Kerr rotation experiments in 3D and 2D n-doped GaAgiSiii and in a p-n junction 
light emitting diode (LED) realized in a two-dimensional hole system^^. However, there is still a debate on the origin 
(extrinsic or intrinsic) of the experimentally observed spin accumulation.i^iHili':^®'^^ 

In this paper we focus exclusively on the theory of the extrinsic effect in a two-dimensional electron gas (2DEG). 
It has long been realized that the extrinsic spin Hall current is the sum of two contributionsi^. The first contribution 
(commonly known as "skew-scattering" mechanism^^'^*^) arises from the asymmetry of the electron-impurity scattering 
in the presence of spin-orbit interactionaSi: electrons that are drifting in the -fx direction under the action of an electric 
field are more likely to be scattered to the left than to the right if, say, their spin is up, while the reverse is true if 
their spin is down. This generates a net z-spin current in the y direction. The second contribution (the so-called 
"side-jump" mechanismi2i22i2i24) is caused by the anomalous relationship between the physical and the canonical 
position operator (see Eq. (|SJ) below). This again leads to a finite spin current in the y direction. The skew-scattering 
and side-jump contributions were widely discussed in the context of the anomalous Hall effecfei24i2iS2422iSi2ii2Si2& in 
magnetic materials. The skew-scattering contribution was first studied by Smi1ii2i2£ while the side-jump contribution 
was introduced by BergeriSSiS^. The theory for both effects has been also discussed recently in several excellent papers, 
both for the extrinsioi^iSi^i and the intrinsic casei^^iSiSSiSi 

In this paper in addition to the previously considered skew-scattering and side-jump contributions, we also include 
the Coulomb interaction effects. The main effect of interactions on the spin transport originates from the friction 
between spin-up and spin-down electrons moving with different drift velocities, the so called spin-Goulomb drag (SCD) 
effecfeSiSiM. We show that while the spin-Hall conductivity associated with the side-jump term is independent of 
the strength of electron-electron interactions, the skew scattering part is reduced by the spin-Coulomb drag, so the 
absolute value of the spin Hall conductivity (and hence the spin Hall current) is reduced for experimentally accessible 
parameters. Since the SCD has been predicted (and recently observed^^) to be a rather significant contribution to the 
overall resistivity in high mobility electronic systems, we think it is important to include it in the description of the 
spin Hall effect, and we show here how this is done. Moreover, we predict that SCD in paramagnetic materials will 
reduce the spin accumulations through the reduction of the skew-scattering resistivity as well as the spin diffusion 
length. Also, we present in the Appendix a simple model for electron-impurity scattering which can be solved exactly, 
leading to an analytical determination of scattering rates as well as side-jump and skew-scattering contributions to 
the spin Hall conductivity. 
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This paper is organized as follows: in Section II the Hamiltonian and the Boltzmann equation are presented; 
in Section III the skew-scattering contribution to the conductivity is derived; in Section IV we use a force-balance 
argument to calculate the side-jump contribution; in Section V the contributions of spin Coulomb drag and spin- 
flip processes are included; in Section VI the spin accumulation in the presence of skew-scattering, side-jump, and 
electron-electron interactions is calculated. We summarize the paper in Section VII. 



II. HAMILTONIAN AND BOLTZMANN EQUATION 

We consider a strictly two-dimensional electron gas (2DEG) that lies in the x — y plane. The hamiltonian is 



H — Hq + Hso + He + He 



where 



2m 



(1) 



(2) 



is the non- interacting hamiltonian (m being the effective mass of the conduction band) , including the electric potential 
Vei{r) generated by randomly distributed impurities. 



2 ^ Cfclr-i 



is the electron-electron interaction (screened by the background dielectric constant e;,), 



Hso ^a^i^p.x \v^Ve^{n) + V iV^^ | • a, , 



(3) 



(4) 



is the spin-orbit interaction (SOI) induced by the electric potential of the impurities, Vei{ri), and of the other electrons 



E 



3=^i i^b\ri-rj \ ' 



and finally 



(5) 



is the interaction with the external electric field E. 

The various spin-orbit terms appearing in the hamiltonian can all be shown to arise from a single basic fact, namely, 
the change in form of the physical position operator under a transformation that eliminates the coupling between the 
conduction band in which the electrons of interest reside, and the spin-orbit-split valence band. If we denote by 
the canonical position operator of the i-th electron, then the physical position operator is given by 



rphysS = "Ti - a{pi X Gi) 



(6) 



and correspondingly the velocity operator is: 



^ \Tphys , H] 



m 



V^Ve^in)+V,V:,+eE 



(7) 



The spin-orbit "coupling constant" a takes into account the effective SOI induced by the valence bands (heavy 
holes, light holes, and split-off band) on conduction electrons in the framework of the 8-band Kane model. Within 
this model one finds^ 



1 



1 



El {E, + Asor 



where Eg is the gap energy between conduction and heavy /light holes bands, Aso is the splitting energy between 
heavy/light holes and split-off bands, P is the matrix element of the momentum operator between the conduction 
and the valence-band edges, and rrie is the bare electron mass. Using values of the parameters appropriate for the 



3 



2DEG in Alo.iGao.gA&ii with a band mass m = 0.074me we find: ah = AAA^ . In this paper we treat the spin-orbit 
interaction to the first order in a, which is justified by the smaUness of the parameter ah/a^, where as — 100 A 
is the effective Bohr radius. Also, we consider the first order corrections from the electric potential associated with 
electron-electron interactions to the spin-orbit hamiltonian, i.e. the terms {ah/a^){e'^ /ebhvp). 

Notice that the canonical positions fi and the canonical momenta pi of the particles are vectors in the x ~ y plane, 
and so is the V operator. Therefore px iVeiin) iV^^)] is a vector in the z-direction, and the spin-orbit interaction 
conserves the z-component of the spin of each electron. This nice feature of our strictly 2D model allows a particularly 
simple analysis of the spin Hall effect without sacrificing any essential features of the spin-orbit interaction. Processes 
that flip the z-component of the spin will be considered separately (see Section V). 

We begin to exploit the conservation of <7z by defining the quasi-classical one-particle distribution function fair, k,t), 
i.e. the probability of finding an electron with z-component of the spin Sz = ^a, with a = ±1, at position r with 

momentum p = hk ai the time t. In this paper we focus on spatially homogeneous steady-state situations, in which 
fa does not depend on r and t (for a discussion of non- homogeneous spin accumulation effects see Section VI). We 
write 

/,(f,fc,t) = /o.(6fc) + /i,(fc) , (8) 

where foa{^k) is the equilibrium distribution function - a function of the free particle energy ek = -^^^ - and fia{k) 
is a small deviation from equilibrium induced by the application of steady electric fields Ea (cr = ±1) which couple 
independently to each of the two spin components. Then to first order in Ea the Boltzmann equation takes the form 

^ hh ■ ^ 

-eEa-—fU^k)^flaik)c, (9) 
TO 

where fia{k)c is the first-order in Ea part of the coUisional time derivative fa{k)c due to different scattering processes 
such as electron-impurity scattering, electron-electron scattering, and spin-flip scattering. As usual, fa{k)c is written 
as the difference of an in-scattering and and an out-scattering terms. For example, in the case of spin-conserving 
electron- impurity scattering one has: 

fa{kUmp = - E [^kk'afS) ~ W~^r^a f ^(k')] 5{eka ~ , (10) 

k' 

where W^p^ is the scattering rate for a spin-cr electron to go from k to /c', and l^a is the particle energy, including 
the additional spin-orbit interaction due the electric field Ea- 

The last point is absolutely vital for a correct accounting of the "side-jump" contribution. We must use 

Efccr = Efc + 2eaha{Ea y- z) ■ k , (11) 



where the second term on the right hand side differs by a factor 2 from what one would surmise from the intuitive 
expression -I- eE • Vphys ■ Why? The reason is that the (5-function in Eq. H10|l expresses the conservation of energy in 
a scattering process. This is a time-dependent process: therefore the correct expression for the change in position of 
the electron /^Tphys must be calculated as the integral of the velocity over time. From the commutator of the physical 
position operator with the Hamiltonian we easily find 



V,V,,{r,) + V,Vi^+eE 



X z 



(12) 



where the term in the square brackets is (minus) the total force acting on the i-th electron. The time integral of this 
term over the duration of the collision (be it an electron-impurity or an electron-electron collision) gives the change 
in momentum Ap during the collision. Thus we see that the change in position is lS.rphys — ^^aa^p x z (this is the 
so-called "side-jump"): hence the change in energy is correctly given by Eq. Hll() . 

Kohn and Luttinger— , have shown that the above form of the collision integral is correct up to third order in 
the strength of the electron-impurity approximation: this is one order higher than the Born approximation and 
should therefore be sufficient to capture the skewdness of the scattering probability, which arises from terms beyond 
the Born approximation. Notice that the collision integral does not contain the tempting but ultimately incorrect 
"Pauli-blocking" factors 1 — faik')- 
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Similarly, the electron-electron contribution to the coUisional derivative has the form24 
Mk)c,e-e ^ - E Wc{ka,p- a; k'aj - a) {/<,(fc)/_,(p)[l - Mk')][l - 

k'pp' 

-Uk')f^,{f)[\ - fAkm - f-.{p)]} Sup.k,+f,S{7k. -?fcv -^p'-a) , (13) 



where Wc (ka, p~ a; k'a, jf a) is the electron-electron scattering rate from ka, p — a to fc'a, — a, and the Pauli factors 
fcr{k), 1 — fa{k') etc. ensure that the initial states are occupied and the final states empty as required by PauH's 
exclusion principle. Notice that, for our purposes, only collisions between electrons of opposite spins are relevant, 
since collision between same-spin electrons conserve the total momentum of each spin component. Accordingly, only 
the former have been retained in Eq. p3fl . 

III. SKEW-SCATTERING 

Let us, at first, neglect the electron-electron interaction. From the general theory developed, for instance, in Refi2&, 
one can easily deduce that the scattering amplitude from one impurity in two dimensions has the form 

4fc',.=^fc£'+^5£fc'(^x^%' (14) 

where A^p and Bj:j:, are complex scattering amplitudes, and k and k' are the unit vectors in the directions of k and 

k' respectively. The second term on the right-hand side of Eq. H14II is due to the spin-orbit interaction. Squaring the 
scattering amplitude and multiplying by the number Ni of independent scattering centers we arrive at the following 
expression for the scattering rate from k to k': 



where 

wh' = M\Ar,'?^\Br,'n - (16) 

and 

WS-, = 2N.,Re Wk'Bl,] . (17) 



Here and in the following = -^^^ is the free particle energy. The second term in the square brackets of Eq. 
depends on the spin (cr = or —1 for up- or down-spins respectively) and on the chirality of the scattering (i.e., 
the sign of (fc x k')z). We will refer to this as the skew-scattering term. It should be noted that for a centrally 
symmetric scattering potential - the only case we are going to consider in this paper - WS-^^ and WS-^^ depend only 

on the magnitude of the vectors k and /c', which are equal by energy conservation, and on the angle 9 between them. 
Furthermore, they are both symmetric under interchange of k and k' - the antisymmetry of the skew-scattering being 
explicitly brought in by the factor (fc x k')z = sin0. Thus, in the following, we will often write Wt^J^ = W^/°'{k,9), 

where W'l''{k,9) are even functions of 9. Notice that the skew-scattering term vanishes when the scattering is treated 
in the second-order Born approximation^^!. Indeed, within this approximation A^p is purely real and B^p is purely 
imaginary, so WS-^, is zero. 

The linearized Boltzmann equation can be solved exactly under the assumption that Wi-.^ and M^^r, depend only 
on the energy = e^i and on the angle 9 between k and k' . The solution has the form 

fa{k) = ha{ek) - fU^k)hk ■ Mk) , (18) 

where Va-{k) is proportional to the electric field. In view of Eq. (|11() it is convenient to expand 

foa{£k) = foaieka) - /q^ (cfe ) (e/ca " Cfe) ■ (19) 
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so that our Ansatz H18(l takes the form 

Uik) = foaieka) - 2f[,^{ek)eaha{E„ x z) ■ k - f[,„{ek)hk ■ V^k) . (20) 

The advantage of this form is that the "zero-order term" foai^ktr) makes no contribution to the coUision integral 
(|10|l . Then, making use of Eq. (|15|l and discarding terms proportional to aW"^ (which are small since W"^ itself is 
proportional to a) we arrive at the following form for the linearized collision integral: 

hAk),^,mp = -^W^|£,{/i.(fc)-/i.(fc')}'^(efe-efe') 

k' 

- "ll^U^^ ^ ^')- {/i-(^) + h'^ik')]5{ek - ek') 

k' 

+ 2aY,WS-f;,^{ek)eah{E, x z) ■ {k - k')5{ek - CfeO , (21) 
fc' 

where fia(k) = — fQ^{€k)hk ■ Vcr{k) is the deviation of the distribution function from unpertubed equilibrium. 

At low temperature, it is not necessary to take into account the full A:-dependence of V„{k), since the derivative of 
the Fermi distribution fQ^{ek) restricts the range of k to the vicinity of the Fermi wave vectors kpa- Thus we replace 
the function Va-{k) by a constant V^, and determine Va from the consistency condition 



Ehk 
— 
TO 

fc 



- hk 

■ — 

m 



hk -• 

/oa(efc)-E— • (22) 



TO 

fc 



Substituting the collision integral from Eq. H21|l on the right-hand side of this equation and moving its last term to 
the left hand side we arrive at 

-2eaY,— \{k~k')WI^^,a{E, x z)f,^{ek)5{ek - efcOl = ^ ^lfc~— ' V^fMSi^k - e^O 

+ 'yJ2^lk'^Ck^k%{k + k')-V,fMKek-ek') , (23) 
fcfe' 

where Ua- = is the density of cr-spin carriers and kp^ is the corresponding Fermi wave vector. The first term on 
the right-hand side of Eq. (|23|l is parallel to Va-, while the second term is orthogonal to it. Then a simple calculation 
leads to the following expression for Va in terms of the electric field: 

e - 2eaa{Ea x z) V„ K x z 
Ea = her — , (24) 

TO Trr Ta T^" 



and its inverse to first order in a is 

- 2ea(j{E„ x z) , (25) 



TO 

where 



Ea - cr-^E„ X z 



1 mA f°° f^"" 

' deefUe) dd W%e,9){l - cos9) , (26) 



(27) 



To- An^ffeFa Jo 
and 

1 mA '■^'^ 



/ dee/^,(e) / d0 ^{€,9)8^^ . 
Jo Jo 



In the above equations epa = 2m Fermi energy for spin a. In the limit of zero temperature the derivative of 

the Fermi function reduces to fQ„{e) — — (5(e — ep^) and the above formulas simplify as follows: 



1 T^o mA '■^"^ 



d9 W''{kp,9){l ~ cos 9) , (28) 
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and 

1 T^o mA '■^'^ 



de W''{kF,0)siii'e (29) 




FIG. 1: Antisymmetric scattering rate Wii^t ~ Jq^ VK"(e, sin^ S in units of riih/m^A (see Eq. 17911 ) as a function of fc^ 
for a model circular well attractive potential Vq = — 5meV and radius a = 9.45nm (described in the Appendix). We choose the 
parameters typical for the experimental 2DEG confined in Alo.1Gao.9As quantum well i.e. density of electrons and impurities 
n2D — rii = 2.0 X lO^^cm"^, m — 0.074me, and mobility ft = O.lm^/Vs. The effective spin-orbit coupling ah — 4AA in 
accordance with—. 



Fig.l shows the antisymmetric scattering rate Jq^ d9 W^"(e, 6*) sin^ 6* calculated numerically using Eq. (|79|l for a 
model impurity potential Eq. (|58|l presented in Appendix, and for the typical experimental parametergy . 



IV. SPIN-HALL CURRENT AND SIDE-JUMP 

The quantity Va obtained in Eq. H25|) . determines the non-equilibrium distribution, according to Eq. (|18|l . We 
now use this distribution to calculate the current density. In order to do this correctly, however, we must remember 
that the spin-orbit interaction alters the relation between the velocity and the canonical momentum. The correct 
expression for the velocity is given in Eq. Q and in the absence of electron-electron interactions takes the form: 

v,^- + 2a[^Ve^ir^) + eE] x a, . (30) 
m 

The second term on the right hand side of this equation contains the net force resulting from the combined action of 
the external electric field and the impurity potential on the i-th electron. This force must vanish when averaged in a 
non-equilibrium steady-state ensemble, since the average value of the momentum must be stationary. Also, the same 
is true in the presence of electron-electron interactions, provided the Coulomb force is included. Thus, we arrive at 
the simple result 

ja = -en^Va- ■ (31) 
Combining this with Eq. (|25|l we obtain the complete relation between electric field and current density: 

-I- 2e^aan„E„ x z . (32) 

The last term on the right-hand side of this expression is known in the literature as the "side-jump" contribution to 
the current density^- . It comes from the use of rather than e^- in the (5- function of conservation of energy - see 



E„ - a E„ X z 
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discussion in the paragraph following Eq. Inverting Eq. (|32|l we obtain the formula for the electric field in terms 

of the current densities: 



(33) 



where 



is the Drude resistivity, 



A. 



Hence the resistivity tensor, written in the 



basis xi, Xj, yi, has the following form: 



P = 



At 



-pf + Atpf pf 

pf -pf + Aipf 

Pf-Aipf pf 



(34) 



The diagonal part of the resistivity reduces to the Drude formula p^ 



D 



as expected. The spin-orbit interaction is 



entirely responsible for the appearance of an off-diagonal (transverse) resistivity. The latter consists of two competing 
terms associated with side-jump (Ao-p^) and skew-scattering (p**), as seen in Eq. Hence our expression for the 

transverse resistivity is different from the expression presented in previous papers (see for example Refi^) where only 
the side-jump contribution appears in the final formulas. Notice that, at this level of approximation, the resistivity is 
diagonal in the spin indices. 



V. SPIN COULOMB DRAG 



Up to this point we have ignored electron-electron scattering processes, as well as scattering processes that might 
flip the spin of the electrons. As discussed in Ref. (^'^) these processes are important because they couple the up- and 
down-spin components of the current density, giving rise to ofF-diagonal elements, p'^f , of the resistivity tensor. The 
Coulomb interaction, in particular, leads to the phenomenon of the spin Coulomb drag (SCD) - a form of friction 
caused by the relative drift motion of spin up and spin down electrons, and the consequent transfer of momentum 
between them^Si^iSSi^. Both Coulomb and spin-flip scattering can be included in our formulation as additional 
contributions to the coUisional derivative fia{k)c (see Ref. ('^^)). In particular, the Coulomb contribution is given 
by Eq. (|13|1 . Substituting Eq. (|20|l into Eq. I|13() and including the first order corrections from electron-electron 
interactions to the spin-orbit Hamiltonian, we arrive at the following expression for the Coulomb collision integral: 

fa{k)c,e-e ^ ' V Wcika.p- a'k'aj - n)[hV„ - hV_, + 2eaha{E, + E^,) x z]{k~k') (35) 

k'pp' 

/o(T(efc)/o-(T(ej9)/o(T(— efc')/o-cr(— ep')'^fe+p,fc'+p''5(^fccr + Cp-o- - f-k'a - ^p'-cr) , 
where T is the temperature, ks is the Boltzmann constant, and we have made use of the identity foai^k) fo-a-{£p)[^ ~ 

focr{ek')][l ~ fo-aif^p')] = [1 - /o<T(efe)][l - fo-a{ep)]foaiek')fo-a{ep') for eka + (p-a - £k'a - ^p'-a = 0. 

Eq. (|35|l is now inserted into the "consistency condition" (|22|l . and the resulting sum over momenta is expressed 
in terms of the spin drag coefficient 7, i.e. the rate of momentum transfer between up- and down-spin electrons, 
according to the formulai24 

n sr-^ ^ ^ (k — k'Y 

7 = Wc{ka,p- a]k'a,]^ - a)- — —-fQa{tk)h-a{^p)foa{-^k')k-a{-^p') 

kk'pp' 

^ '^£+j?,fe'+p'^(efe^ + ep-<T-efc'<7-ep'-o-)- (36) 
Then equation for Va is modified as follows: 

Ecr - 2 = h cr Jo- H j^a- + 27 eaa(E^ + E^„) x z + — (37) 

m Ta Ta T§ Ua-ne ne n T^ria-e 

The third and fourth terms on the right hand side (r.h.s.) of this equation are connected with the SCD and 7 is 
the spin-drag coefficient calculated in Refi^. The fifth term on r.h.s. of Eo. 1371 comes from the side-jump effect in 
Coulomb scattering. The last term in Eq. H37|) is associated with spin-flip collision processes, characterized by the 
spin relaxation time r^. We include it for completeness, even though spin-flip effects are expected to be small in 
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n-doped semiconductors. The current is now given by the full velocity operator of Eq. Q, but thanks to the force 
balance condition in the steady state it reduces again to the simple form 



Jo 



-en„Vrr 



(38) 



v 



Combining fEa. 137(1 with fEo. 138(1 leads to the resistivity tensor, which in the basis of a:|, xi, m has the form 

p 

-fj - Pi ^ P ti.yti.1 -Pi ^' ^ 



pT' 



^p"" - p\ 



\ 



HP I 



AT 



(39) 



where p^^ — m'^/ne'^ is the spin Coulomb drag resistivity and p'^ = m/riae'^T'^ (recall that = ^^r^ is a dimensionless 
quantity). AJ" and -BJ" represent the terms of the first order in electron-electron coupling 7 and in SO coupling 
a and are defined as follows: A'^" = —\^psDn-a/n„ + 2m*a'y[—n-crP^ /n + (n-ajn — r?_^lnn„)p^^\ and = 
XaPSD + 2m*a'~i[—n-aP-crl''^ + (ri-ajn — Ur, ln)p^^\. Notice that the resistivity satisfies the following symmetry 
relations: 



(40) 



pT. - P'-L (41) 

where upper indices /3 and /?' denote directions, and the lower ones spin orientations. From Eq. ((39(1 one clearly sees 
that the spin Coulomb drag and spin-flip processes couple the spin components of the longitudinal resistivity. Further, 

the spin-Coulomb drag corrections to the spin-orbit hamiltonian (07 corrections) result in the p^„^_„ terms, i.e. in 
the transverse resistivity which couples opposite spins. Also, the SCD renormalizes the longitudinal spin-diagonal 
components p^^„^_„ of the resistivity, in such a way as to satisfy Galilean invariance. 

The spin-flip collision processes come as a phenomenological term which could have origin for example in some 
random magnetic field, which docs not appear in the original spin-conserving Hamiltonian (Eq. (^l). The spin-flip 
relaxation time, is given by the microscopic expression^ 



ddW'J_^{ek.O)cosB . (42) 



where W^'^^ denotes the spin-flip scattering rate from spin a to the opposite spin orientation Since the 



Tl 

relaxation time for spin-flip processes r' is very long,"^^ the SCD normally controls the coupling between the spin 
components, except at the very lowest temperatures (the spin drag rate 7 vanishes as In T while the spin-flip rate 
remains constant). Spin-flip processes will therefore be omitted henceforth. 

In a paramagnetic material there is symmetry between up and down spin densities, mobilities, etc. and one can 
easily separate the spin and charge degrees of freedom. Then combining the Eq. ((37(1 and Eq. ((38(1 simplifies to: 

Ec = p'^Jc + 2ip'' - Ap^ - Xpsd)1 X z , (43) 



Es = 4(p^^ + p^)j, + 2{p'' - Ap^ - XpsD)jc X z 



(44) 



2 ■ 



where we omit spin-flip processes, and the charge/spin components of the electric field are defined as Ec = 
Es — E^ — El, and the charge and spin currents are jc = j| +ji and js = "'^ 2''^ ' respectively. The spin-Coulomb drag 



renormalizes the longitudinal resistivity only in the spin channel. This is a consequence of the fact that the net force 
exerted by spin-up electrons on spin-down electrons is proportional to the difference of their drift velocities, i.e. to the 
spin current. Additionally, the electron-electron corrections to the spin-orbit interactions renormalize the transverse 
resistivity in the charge and spin channels, so the Onsager relations between spin and charge channels hold. Under 
the assumption that the electric field is in the x direction has the same value for spin up and spin down we see that 
Eqs. ((43(1 and l(44() yield the following formula for the spin current in y direction: 



p-/2(p 



D\2 



l + pSD/pD 



X ' 



Ex 



(45) 
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FIG. 2: Spin Hall conductivity as a function of temperature. cTy^open squares), ay% (open and close circles) are the side-jump 
and skew-scattering contributions, respectively, to the conductivity in the absence of electron-electron interactions, al°J is the 
total spin conductivity when electron-electron interactions are taken into account. We choose the parameters typical for the 
experimental 2DEG confined in Alo.1Gao.9As quantum well i.e. density of electrons and impurities n2D = rii — 2.0 x lO^^cm"^, 
m — 0.074me, and two sets of mobilities and relaxations times: ft — O.lm^/Vs, r = 4 x 10~^ns Tbb = 0.02ns and p, = Im^/Vs 
r = 4 X 10~*ns Tbb = 0.2ns. The effective spin-orbit coupling ah — AAA in accordance with—. We used the model potential 
(see appendix) where an effective impurity radius a — 9.45nm, the height of attractive impurity potential Vq = — 5meV for 
fj. = O.lmVVs and Vo = -1.6meV for p. = ImVVs. 



The first term in the square brackets is associated with the skew-scattering, while the second is the side-jump contri- 
bution. Notice that the side-jump conductivity — —Oine^ depends neither on the strength of disorder nor on 

the strength of the electron-electron interaction. In contrast, the spin-Coulomb drag reduces by factor 1 + p^^ j 
the skew-scattering term, so the total spin Hall conductivity and the spin Hall current are reduced as well. 

The temperature dependence of the spin conductivities Uyx = j'^y/E^ for two different mobilities and for the 
parameters of the recent experiment on the 2DEG in Alo.iGao.gAs^^ is presented in Fig. 2. We used Eqs. (|78(l and 
(|79|l in the appendix to calculate the scattering rates. The procedure is as follows. First we find the scattering 
rate t^{= t-„) from the mobility. Then using Eqs. H78I) and (|26|1 we estimate the strength of impurity potential Vq- 
Finally, we compute the skew-scattering rate by numerically evaluating Eq. (|79|) . In general the skew-scattering rate of 
Eq. H26I) is temperature dependent through the Fermi distribution. However, in a wide range of mobilities the impurity 
potential Vq is much smaller than the Fermi energy and the scattering rate does not change with energy around Ep^ 
i.e. [dW°'{€,Q) / de]^=EF ~ (see Fig.l). Since the temperature dependence of the skew-scattering rate comes from 
the energy dependence of 14^° (e, 8), we see that the temperature dependence of skew-scattering term is very weak. 
In Fig. 2, cr^-j, is the side-jump contribution to the conductivity in the absence of electron-electron interactions (found 
from the last term of Eq. I32() . cry% is the skew-scattering conductivity in the absence of electron-electron interactions 
(found from second term of Eq. (|^ ). (t*°* is the total off-diagonal conductivity modified by the spin Coulomb drag 
term, as given by Eq. I|45|l . For T = the total spin conductivity is the sum of two contributions: skew-scattering 
proportional to pss/'ip\) and side-jump proportional to A/2p£) = —ane^. For an attractive electron-impurity potential 
(Vb < 0) the two contributions have opposite signs, consistent with previous theoretical estimates in Refs,^^'^^. The 
ratio of the skew-scattering to side-jump terms depends on the mobility. The skew-scattering conductivity scales as 
p while the side-jump is independent of mobility. As a consequence the skew-scattering conductivity dominates for 
high mobilities while the side-jump dominates at low mobilities. The different ratios of skew-scattering to side-jump 
conductivities reported in recent theoretical papersiii& result from choosing different mobilities (see also Fig. 2). The 
spin-Coulomb drag is the only temperature dependent contribution in this calculations: it scales as InT for T < Tp 
and reduces the absolute value of spin conductivity and spin current. Moreover, the reduction of spin-Hall effect by 
spin-Coulomb drag depends strongly on the Drude resistivity. Hence the reduction of spin-Hall conductivity will be 
of the order of a few percent for low mobilities (invisible change in Uyx for p = 0.1m~/Vs in Fig. 2) and of the order 
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of 25%-50% for high mobihties (see ayx for fi = Im /Vs in Fig. 2). 

Thus far, as stated in the introduction, we have only considered the extrinsic spin HaU effect. What about the 
intrinsic effect? In recent experiments by Sih et aliAi on a 2DEG confined in an (110) AlGaAs quantum well, three 
different contributions to the spin Hall conductivity are present in principle: the impurities, the linear-in-A; SO Rashba 
field, and the cubic-in-fe SO Dresselhaus field which is perpendicular to the Rashba field. Since in (110) quantum 
wells, the Dresselhaus field is in the [110] direction; we do not expect any spin-Hall current or spin-Hall accumulation 
connected with this term in the plane of the quantum well. Further, it has been established that the spin-Hall 
conductivity in an infinite 2DEG with only linear spin-orbit interactions of Rashba or Dresselhaus type vanishes 
for arbitrarily weak disorderi^Si^i*^2ii^. Hence we believe that the extrinsic contribution ay^ ~ 4 x 10~'^l/k^l is the 
dominant one in the referenced paper. Also, the theoretical estimates of extrinsic and intrinsic contributions in 
the experiments on three-dimensional n-doped GaAs— show that the extrinsic contribution is an order of magnitude 
larger than the intrinsic one. 

VI. SPIN ACCUMULATIONS 

Let us now study the influence of spin-Coulomb drag on the spin accumulation. This is relevant to the interpretation 
of the experiments by Sih et ali^ and Kato et ali^. We consider a very long conductor in the form of a bar of length 
L in the x direction and narrow width W in the y direction. A charge current flows only in the x direction. The 
y-components of the current j-^, with a =| and | add up to zero everywhere and individually vanish on the edges 
of the system, i.e. = at y — ±W/2. In order to satisfy the boundary conditions the system cannot remain 
homogeneous in the y-direction. A position-dependent spin density, known as spin accumulation develops across the 
bar, and is reflected in non- uniform chemical potentials fia{y)- In the steady state regime the spatial derivative of 
the spin-current in the j/-direction must exactly balance the relaxation of the spin density due to spin-flip processes. 
This condition leads to the following equation for the spin chemical potentiak'^'i 

d^y ^ Ll ' 

where is the spin diffusion length. The solution of this equation is: 

M.(y) - fi^Ay) = Ce^ -f C'e-i , (47) 

where C, C' are constants to be determined by the boundary conditions j\„{±.W /2) = 0. The effective electric field 
in the y- direction is eE^ = d^a/dy. Thus, 

-e{Ey - El) = ^ ^e-v/i. _ ^^v/l. (48) 

dy Lg Lg 

Using the boundary conditions at y = ±W/2 we obtain 

Eyi±w/2)^py^^j^ + py^^jf m 
Efi±w/2) = py^,j^ + pyijt (50) 

Making use of Eqs. (|48I50|) to eliminate the electric field, we obtain the following set of equations: 

f e^^/^^= - f = f - f g-^/^^' = -e[pf^jf + pf J=jf - p^j^ - pUf,] , (51) 

Ijg L/g L/g L/ g 

which gives immediately the solution C ~ ~C' and 

2cosh{W/2Lg) ■ ^ ' 

Thus, the formula for the spin accumulation is 

eLg[pY^j^ + pIIjI - pf^^jf - plljf] smh{y/Lg) 
MT (y) - Mi (2/) = cosh(M^/2L,) ' ^^^^ 
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Finally, upon substituting the matrix elements of the resistivity from Eq. (|39|l we find 



jf (pf^ - Atpf + + B^") + jlipf - Axpf + A]^ + BJ") sinh(y/L,) 



- = cosh(M^/2L.) • 

The spin-Coulomb drag modifies the spin accumulation in three different ways: (i) the spin Coulomb drag resistivity 
appears directly in the terms A'^°' and B^", defined after Eq. (jSHJ; (ii) the values of the spin components of the 
longitudinal current are, in general, affected by the SCD; (iii) the spin diffusion length Ls is significantly reduced 
by the spin Coulomb dragS^ as seen from the formulaS^. 

Ls = , (55) 

1 + PsdIpd 

where xo is the susceptibility, xs is the spin susceptibility, and Lc is the density diffusion length. 
In paramagnetic materials Eq. 154|) simplifies to: 

. ^ , . _ -2eLsjx[p"' - - Xpso] sinh(y/is) _ -2eLsEx[p"' ~ Xp^ - Xpso] smh{y/Ls) . 
p^(y)-p,{y)- cosh(W^/2L,) " pz5 cosh(P^/2i,) ' ^^^> 

where jx — E^/pd indeed is independent of psD- Further the Eq. (|56|l at the edges of sample for L — W/2 gives: 

piiiW/2L,) - Mi(T^/2L,) - -2eL,2x¥' ~ Ap^ - Xpsu] tanh(iy/2L,) . (57) 

The three terms in the square bracket of Eq. H57|l are the skew-scattering term, the side-jump contribution and the 
electron-electron correction. The last term reduces the spin accumulations. Additionally, the spin-Coulomb drag 
reduces the spin accumulation through the spin diffusion length (see Eq. However, in the limit oiW ^ Ls^ 

tanh(W/2Ls) can be approximated by W l2Ls^ and the spin accumulation at the edges becomes independent of Lj. 
In this limit, the influence of SCD on the spin accumulation is only through the XpsD term. Notice that in the limit 
of infinite spin-relaxation time {Ls oo) the spin accumulation can be obtained directly from the homogeneous 
formulas, Eqs. (gSJ and 

For a two-dimensional electron gas confined in Alo.iGao.gAs quantum well measured by Sih et ali^ with electron and 
impurity concentrations = 7120 = 2 x 10^'^cm^^, mobility /2=0.1m^/Vs, Ls = 1/im, t — i x lO^^ns, Tss — 0.02ns, 
ah = AAA and for the sample with width W = 100/im, we calculate the spin accumulation to be — 1.5meV/|e| on 
the left edge of the sample (relative to the direction of the electric field) i.e. for W/2 — —50pm. This means that 
the non-equilibrium spin-density points down on the left edge of the sample and up on the right edge. It is not clear 
at present; whether, this sign of the spin accumulation is consistent or not with the sign of the Kerr rotation angle 
measured in the experimentsiS*ii. 



VII. SUMMARY 



We have developed the microscopic theory of the extrinsic spin Hall effect taking into account the skew-scattering, 
side-jump and Coulomb interaction effects. The total spin conductivity in zero temperature is a sum of the skew- 
scattering and side-jump terms. The spin Coulomb drag is the only temperature dependent term, causing a quadratic- 
in-T reduction of the spin Hall conductivity. Further, we find that the spin Hall conductivity associated with the 
side-jump contribution is independent of the strength of electron-electron interactions, while the part of the spin 
conductivity connected with the skew scattering is reduced by the spin-Coulomb drag for experimentally accessible 
mobilities. Moreover, we predict that in paramagnetic systems the spin-Coulomb drag reduces the spin accumulations 
in two different ways: (i) directly through the reduction of the skew-scattering contribution (ii) indirectly through the 
reduction of the spin diffusion length. 
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IX. APPENDIX- AN EXACTLY SOLVABLE MODEL FOR SKEW-SCATTERING 



We present the calculations of the scattering rates for a circular well potential of the form: 

V{r) = Vae{a - r) + aaL,S,S{r - a)Vo , (58) 

which is attractive for Vb < and repulsive for Vq > 0. The second term on the right is the spin-orbit interaction 
and a = ahja? where a is an effective spin-orbit coupling found in 8-band Kane model and a is the impurity radius. 
Since the orbital angular momentum Lz = I and the spin angular momentum Sz = <J are conserved we can separate 
the wave function into radial and orbital parts: 

^kUr,e)^Rua{r)e'"^ (59) 
and the corresponding Schrodinger equation has the form: 

- ^ (r'L + -R'kia] + P-^Rkia + V{r)Rki. (r) = ERki. (r) (60) 

where E = h^k"^ I2m and the prime denotes a derivative with respect to r. We now express lengths and wave vectors 
in units a and respectively, so r should be understood as r/a, k as fca, and, of course, a = 1 in these units. The 
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dimensionless Schrodinger equation is 

R'L + \R'kia + f - ^0 - 5 ) Rki. = 0, r < 1 (61) 

and 



1„, P 

r.2 



r \ r 



Rli, + -R',i, + fc' - ^ Rkia - 0, r > 1 , (62) 



where vq = ^™^»° ig the dimensionless parameter which measures the height of the impurity potential barrier. The 
regular solution of this equation for r < 1 is 

R.kia{r) = Jmii^r), r<l, (63) 

where we have defined i' — \/k'^ — vq. On the other hand, the solution for r > 1 can be written as a superposition of 
the two independent solutions of the differential equation H62|) : 

Rkia{r) = e'^''[cos5i^J\i\{kr) - sin5i„Y\i\{kr)],r > 1 (64) 

The matching conditions on the wave function and its derivative lead to the following equations: 

Rkia{l+) = i?fci.(l") = RkiM) (65) 



KiA^+) ~ KiA^-) = alaRkia{l)vo (66) 
Substitution of Eqs. (|63|l and H64|) to the matching conditions yields 

k cos 5i„J[i\{k) - k sin Sia-Y^^{k) J\i\{'^) 



cosSiaJ\i\{k) - sin(5/<jY|/|(fc) J\i\ii^) 
from which one gets the following equation for the phase shifts Sia-: 



alavo , (67) 



fci^;i(fc)-A.i|^i(fc) 

where /3ia- — v + alavo, J|;|(fc) and Y^\{k) are the Bessel functions of the first and second kind. The wave 

function at large distance from the scattering center can be written as: 

i^kiAr, 0) a/^cos (fcr - - 5 + 5,<.)e^'^e^^- = ^kiAr, 0) + £!:^e^('c'-l'IV2-V4),i^« (69) 
V Trkr 2 4 V27rfcr 

where ip^i^ir, 0) — ^ cos (fcr — — j)e*'^ is the free wave function in the channel of angular momentum /. The 
scattering amplitude fa{k,6) is the factor multiplying the outgoing wave in the above equation: 

c=o 2iSi 1 
/ — — oo 

The differential cross section is accordingly given by: 

= 0)|2 = _L V(e2'''" - l)(e-2'5,'. _ i)e— /2(|i|-|r|)g»(i-r)e (7^) 

fit/ / _ An r\! 

' ° 1,1' 

Finally we notice that the total scattering rate is related to the differential scattering cross section for a single 
impurity as follows: 

W{k, 9) = W'{k, 9) + aW%k, 6) sm9 = , (72) 
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where = Ni/A is the areal density of impurities. Combining this with Eq. (|71|l we find 

2^2 ft3 



W{k,e) = n. 



(73) 



To identify W and we separate Eq. (|73ll into even and odd components with respect to the scattering angle 0, 
which can be easily done using the identity e**^'^' — cos[(Z — l')d] + i sm[{l — l')9]. Then 



271^ 

nrfiA 



^(g2»5,„ _ i)(g-2.Av„ _ cos[(/ - i')e] 



LI' 



and 



m^A sin 9 



LV 



Making use of the identities e^^'*'" — 1 = — .^'^V- we rewrite the scattering rates as 

^ ' ' * mM ^ (cot(5,^ - i)(cot(5;'^ +i) 



and 



W''{k,9) = crn, 



■E 



TO^^ (cot (5;cr — i) (cot Si'c, + i) sin ( 



(74) 



(75) 



(76) 



(77) 



where the phase shifts are completely determined by Eq. (|68|l . Notice that the sums over I and I' in Eq. (|75|l and 
Eq. (|77() run from —oo to oo and the phase shifts have the symmetries 5-i^-a(oi) = 5i_cr{oi) and 5_;.(j(— a) = 5i_cr{ot), 
which implies that W^{k^ 9) and W°'{k, 9) are invariant under spin reversal a — > —a and M^°(fc, 9) changes sign with 
a change of sign of a, as expected. The integral over 9 (Eq. I28|) eliminates the majority of the terms from the sum 
over I and V and the only non-zero terms for W^{k^ 9) are with I — I' and I = I' ±1. This gives 



de W'{k,9){l-cos9)^n,^—r 

m'^A 



E 



E 



^ (C0t(5f^ + 1) ^Jp^^ (C0t(5;^ - i)(cot(5;'CT + «) 



(78) 



One can see that W^''(fc, 0) is modified by spin-orbit interactions, however it has a non-zero value even if spin-orbit 
interactions are absent. We checked numerically that the sum over / in (Eq. I78|) is convergent after taking into 
account a few first terms. For the skew-scattering rate W°'{k,9), the integral over 9 fEa. I29f) have non-zero terms 
only if Z = /' ± 1. This gives: 







2Tr «7r3fi3 

deW''{k,9)sin^9^an,--^ ^ 



■m?A ^^p^^ (cot^io- - i)(cot(5;'cr + «) 



(79) 



For very small ka the only relevant terms are /' = 0, |^| = 1 and \l'\ = 1, Z = which yields to: 



/ d9 Wik, 9) sm^ 9 -an., 
Jo 



87r3ft3 
m'^A 



1 



1 



(cot Sic — i)(cot Sqo- + i) (cot (5ocr — i)(cot + i 
1 1 



1671^^3 

■ m?A 

16712^3 



Re 



(cot^Ocr — z)(cot(5iCT + i) (cot(5_icr — i)(cot(5ocr + i) _ 
1 1 



(cot 5la - i)(cot (5ocr + i) 

1 -f cot (5o(T cot 5i„ 



(cot - i)(cot(5ocr + «) 

1 + cot 5{)„ cot (5_icr 



(cot^ -f l)(cot2 5a^ + 1) (cot^ (5_i^ + l)(cot2 5q„ + 1) 



(80) 



independent of 6i. Notice that ^6* Wilt, 9) sin vanishes in the absence of spin-orbit interactions, since all the 
phase shifts are independent of a in that case. 



